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ON  THE  STUDY  OF  THE  SHELF-LIFE  OF  COMPONENTS  OF  AUTOMATIC  SYSTEMS 


[Following,  is  a  translation  of  an  article  by  G.  V.. 
Druzhinin  in  Izvestiya  Akadetnii  Nauk,  otdelenie  tektinic- 
heskikh  nauk,  Energetika  i  Automat ika  (News  of  the  - 
Academy  of  Sciences,  Department  of  Technical  Sciences, 
Power  Engineering  and  Automation),  No.  6,  1959,  Moscow, 
pages  141-150]  ,,  *1  ’ ‘ 


The  peculiarities  of  studying  problems  of  shelf  life  of  auto¬ 
matic  elements  are  examined  in  this  article.  For. the  study  of  shelf 
life  and  of  the  reliability  of  the  elements,  it  is  proposed  to  uti¬ 
lize  linear  random  functions,  all  realizations  of  which  are. straight 
lines,  :■  The  properties  of  these  random  functions  are  described 
briefly.  ,The  problem  is  considered  concerning  the  computation  of  the 
basic  characteristics  of  the  shelf  life  of  automatic  elements  •.--  ,of 
the  law  of  distribution  of  the  hold-up  time  in  operable  condition, 

A  number  of  the  most  important  properties  of  automatic  ele¬ 
ments  is  related  shelf -life,  that  is,  the  capability  to  be  found  in 
proper  condition  in  the  process  of  storing.  In  other  words,  shelf- 
life  is  the  reliability  of  non-operating  (stored)  elements  and  systems. 
One  can  usually  consider  the  operation  of  elements  (found  under  cur¬ 
rent  or  voltage)  as  a  particular  case  of  storage.  • 

For  the  qualitative  measure  of  shelf-life,  let  us  apply,  by 
analogy  with  the  common  qualitative  determination  of  reliability, 
the  probability  s  of  finding  the  element  or  the  system  in  an  operable 
condition  in  the  course  of  the  assigned  time  t  of  storage  under  ,  ; 
definite  conditions.  We  will  call  the  relationship  s(t)  the  shelf-life 
function.  The  shelf -life  function  s  (t)  or  the  theory  of  .probability 

of  hold-up  time  in  proper  condition  is  f(t)  *  -  ds(t)  or  the  breakdown 

.  1  ....  :  ■  ;.-y  ;  i  ,  ,  ......  _  dt  v.  " 

rate  A  (t)  =  f (t)  characterize  the  shelf  life  of  an  element  or  a 

s(t)  V.v  V: .V  3 

system.  ,  ..  ..  ..  .  ;  .,  .. . 1  . .  .... 

Prediction  of  the  shelf-life  consists  of  computation  of  one 
of  these  characteristics. 

.1.  The  determining  parameter  of  an  element  as  a  random  function 

of  time.  During  storage,  the  basic  bulk  of  failures  occurs  owing  to 

aging  of  the  elements.  Every  element  has  the  determining  parameter); 
which  can  serve  as  a  measure  of  its  quality,'  The  determining  para-  :> 
meter  of  an  element  occasionally  changes  in  the  process  of  storage 
(the  element  ages).  ,  •  . ,  .  .  .  "  "  V. .. 


During  aging,  a  variation  of  the  determining  parameter  of 
an  element  usually  occurs  monotonously  and  in  one  direction  for  all 
elements  of  a  given  type.  For  example,  all.  carbon  composition 
resistors  monotonously  decrease  their  conductivity  with  the  passage 
of  time. 

The  determining  parameter  of  an  element,  during  its  change, 
can  achieve  a  critical  value  during  which  the  condition  of  the  element 
is  considered  unsatisfactory,  that  is,  the  element  goes  out  of 
order  [breaks  down].  This  critical  valuei  of  the  determining  para* 
meter  is  established  as  a  constant  for  a  group  of  similar  dements 
which  are  in  storage  and  which  are  not  connected  to  a  system.  If 
one  considers  the  devices  (systems)  which  consist  of  many  elements, 
then  the  critical  value  of  the  parameter  *1  for  each  element  will 
be  determined  by  the  location  of  the  element  in  the  skeleton  diagram 
of  the  system.  For  brevity's  sake,  the  critical  value  of  the 
determining  parameter  we  will  sometimes  call  the  limit  of  the  element. 

Breakdowns  (failures)  of  elements  owing  to  aging  are  called 
progressive  which  indicates  the  character  of  their  origination.  The 
moment  of  appearance  of  a  progressive  -failure  can  be  forecast.  This 
can  be  achieved  either  by  means  of  observation  of  the  behavior  of 
each  specific  elements  (boundary  control  [l])  or  the  problem  can  be 
examined  statistically.  The  determining  parameter  of  the  element 
can  be  considered  as  a  random  function  of  time  H  (t).  At  a  fixed 
moment  of  time  t*,  the  value  of  the  determining  parameter  will  be 
the  random  value  H^.  In  principle,  H  (t)  can  be  considered  as  a 
vectorial  random  function  of  time.'  The  below  presented  results  are 
easily  generalized  for  the  case  Of  an  n -measured  random  vector. 
However,. in  practice,  there  does  not  usually  arise  the  need  for  such 
a  generalization. 

In  the  process  of  storage,  along  with  gradual  failures,  a 
certain 'number  of  sudden  failures  of  elements  is  observed.  During 
sudden  failures,  there  takes  place  a  sharp,  practically  instantaneous 
changes  of  the  elements 'characteristics.  The  discontinuities  of 
y,  j(t)  of  the  random  function  H(t)  correspond  to  sudden  failures  of  . 
the’* elements.  The  time  of  occurrence  of  the  railization  discontinuity 
y,  ,(t)  is  a  random  value.  Owing  to  this,  the  random  function  H(t) 
is  continuous  even  in  the  presence  of  sudden  failures.  It  follows 
from  the  determination  [2]  of  the  continuity  of  a  random  function. 

The  determining  parameter  of  an  element  can  always  be  chosen  so  that 
it  would  decrease  with  the  passage  of  time  and  during  a  sudden 
failure  ■  0. 

!.'■  In  accordance  with  the  stated  [remarks],  a  sudden  failure  can 
be  considered aa  individual ' occurrence  of  the  change  of  the  elements' 
parameter  and  correspondingly  can  proceed  during  the  computation 
of  the  characteristics  of  the  random  function  H(t)  according  to  the 
experimental  data.  Under  conditions  of  storage,  it  is  practically 
impossible  to  achieve  continuous  control  of  the  elements'  adjustment. 


Information  concerning  sudden  failures  is  obtained  only  periodically 
at  the  moments  of  measurement  of  the  elements'  parameters.  Therefore, 
during  computation  of  the  characteristics  of  the  random  function 
H(t)  according  to  the  experimental  data,  it  ia  necessary  to  consider 
as  equal  to  zero  the  values  of  the  determining  parameters  of  the 
elements  which  suddenly  fail  at  a  given  moment  of  time  in  storage. 

Thus,  in  the  general  case  of  the  presence  of  both  types  of 
failures,  one  can  utilize  methods  developed  under  the  supposition 
that  all  failures  of  elements  are  gradual. 

The  assumption  concerning  the  normal  law  of  distribution 
fi(Vj  )  of  the  values  of  the  determining  parameter  of  an  element  at 
any  fixed  moment  of  time  t*  are  in  good  agreement  with  the  experimen¬ 
tal  data.  Therefore,  further  discussion  will  frequently  be  conducted 
under  the  assumption  that  the  random  function  H(t)  is  distributed 
normally. 

Study  of  the  shelf-liSe  of  elements  by  means  of  the  examina¬ 
tion  of  random  functions  H(t)  has  the  following  features: 

1 .  The  variation  of  the  determining  parameter  of  each  specific 
element  during  the  course  of  storage  time  can  usually  be  approximated 
as  a  linear  function  with  sufficient  accuracy. 

2.  Usually  an  element  ages  morally  [moral 'no]  more  rapidly 
than  physically.  Therefore,  any  study  of  shelf-life  is  always 
associated  with  the  extrapolation  of  the  determining  parameter  of 
the  element.  The  forecast  of  the  shelf-life  is  principally  disting¬ 
uished  by  this  from  a  forecast  of  the  reliability  or  the  accuracy 
of  elements  and  systems  of  transitory  effect,  where  the  results  of 
tests  of  elements  statistically  characterize  their  properties.  Linear 
extrapolation  is  more  simple  and  convenient. 

3.  The  opportunities  of  conducting  measurements  of  the  values 
of  the  elements'  determining  parameters  are  very  limited.  In  many 
cases,  in  the  process  of  storage  or  exploitation,  it  is  possible 

to  measure  only  once  the  values  of  the  determining  parameter  of 
similar  elements.  r  ... 

2.  Concerning  the  extrapolation  of  the  determining  parameter 
for  a  minimum  number  of  measurements.  The  above  indicated  features 
of  the  study  of  the  shelf-life  of  elements  define  the  necessity  of 
applying  linear  extrapolation  of  an  element's  determining  parameter. 

In  addition,  it  is  convenient  to  consider  the  determining  parameter 
H(t)  as  a  linear  random  function  of  time,  that  is,  as  a  random 
function,  all  realizations  of  which  are  straight  lines. 

A  linear  random  function  of  time  H(t)  is  the  linear  function 
of  two  random  arguments:  of  the  initial  ordinate  H0  and  of  the 
angular  coefficient  B  of  a  straight  [line] 

H(t)  =  H0  +  Bt  (2.1) 

The  expected  value  of  a  linear  random  function  is  .  •. ! 

m^  (t)  =  +  mb  t  (2.2) 


where  m>,^  —  is  the  expected  value  .of  the  Initial  ordinate  Hq;  tn^ 

--  is  the*  expected  value  of  the  angular  coefficient  B, 

The  correlation  function  of  a  linear  random  function  is 


Kjt,?)  +  it  +t')«h  +tt% 


(2.3) 


Here,  Dv,  ,  ,  Db  —  is  the  variance  of  the  initial  ordinate  Hq 
and  of  the  angulhr  coefficient  B,  .  —  is  the  coupling  moment 

of  these  random  values.  r  ' 

Assuming  in  formula  (2,3)  that  t  *=  t\  we  will  obtain  the 
formula  for  the  variance  of  a  linear  random  function 

Ori(t)  *  DVc  f  Zihj,ne  i-tOt,  (2.4) 

Utilizing  formula  (2.4)  one  can  represent  formula  (2.3)  as 


KAttv)  *Dnit)  f  £  OkJfl.  [t  ~t') 

d  *- 


(2.5) 


Thus,  the  variance  of  a  linear  random  function  determines  the 
form  of  its  correlation  function.  If  the  linear  random  function  is 
disturbed  normally,  then  the  expected  value  mn  (t)  and  the  variance 
(t)  are  its  comprehensive  characteristic. 

One  can  assign  constant  numerical  characteristic  to  a  normal 
linear  function.  One  can  be  convinced  of  this,  having  examined 
formulas  (2.2)  -  (2.4)  where  all  coefficients  are  constant. 

The  expected  value  and  the  variance  of  a  linear  random  function 
can  be  determined  in  accordance  to  the  statistical  data  of  two  of 
its  sections.  Actually,  in  agreement  to  formulas  (2.2)  -  (2.4),  to 
determine  the  expected  value  ffl>,  (t)  and  the  variance  Dv»  (t)  (or  the 
correlation  function  (t,  t')  of  the  linear  random  function  H(t), 
it  is  necessary  to  know  the  five  numerical  characteristics:  the 
expected  values  mw  and  m^  ,  the  variances  and  Db  and  the 

correlation  moment  >ip  .  The  numerical  characteristics  of  these 
cross  sections  t^  and  t^  +  the  expected  values  m^('  and  mn;f/  , 
the  variances  D  and  D*,.  t  and  the  correlation  moment  of  the 
cross  sections  k„>t\U'» 

Let  us  find  the  relationships  of  the  numerical  characteristics 
of  a  linear  random  function  to  the  numerical  characteristics  of  two 
of  its  sections  t^  and  t^  +  i 


The  angular  coefficient  of  the  straight  line  is 


Hi -  lj_j 


The  initial  ordinate  is 


//, * 


(2.6) 


(2.7) 


l  +  i 


In  accordance  with  formulas  (2.6),  and  (2.7) 


M  id]  s 


to 


V/ 


^  t  +  i  ;■**  1  ( 


hi,,:  -  jV 


(2.8) 

(2.9) 


'^Utilizing  formulas  (2.6)  -  (2.9),  one  can  obtain  expressions 
for  the  remaining  characteristics  of  the  linear  random  function  H(t) 


K 


X  H 


t£>»;»i  h 


(2.10) 


a 


i .*! 


-*i> 


Dh  r  ■  +'On£  ■  (2.11) 

0y}t  2  "  (2.12) 

<x  L*i  -tj  ** 


.  Thus,  to  determine  the  characteristics  of  a  linear  random 
function  of  time  from  the  test,,  it  is  necessary  to  measure  the  values 
which  this  function  takes  at  only  two  instants  of  time.  Moreover, 
since  the  initial  values  of  the  realizations  of  a  random  function 
(consequently  also  the  numerical  characteristics  m*t  and  )  are 
known  in  practice,  then  usually  the  problem  of  determining  the 
characteristics  of  the  linear  random  function  H(t)  results  in  the 
measurement  of  the  values  of  the  random  variable  at  one  of  the 
instants  of  time  t^  and  in  the  computation  according  to  these  results 
of  the  numerical  characteristics  D^y  and  K^.  Having  established 
in  formulas  (2.8),  (2.10)  and  (2.11)  t£  -  0  ana  tj,  +  i  *  ti,  one  can 
compute  according  to  available  data  mjj,  D^,  and  kj,  ^ 

Sometimes,  depending  on  the  conditions  of  the  problem,  measure¬ 
ment  of  the  values  of  the  various  realizations  rt j  (t )  of  the  fandom 
function  H(t)  at  various  times  (Figure  1)  is  found  to  be  expedient. 

If  the  initial  value  of  the  specific  realization  is  known,  then  it 
is  necessary  to  conduct  still  Snother  measurement  at  any  sufficiently 
remote,  from  zero,  moment  of  time.  During  such  a  method  of  measure¬ 
ment,  the  values  and  b  ate  individually  computed  for  each  realiza 
tion.  The  obtained  data  makes  it  possible  to  compute  the  statistical 
numerical  characteristics  m ,  mjj,  D ,  D^,  and  kjj>^  in  accordance 
to  known  formulas  of  the  theory  of  probability. 

Let  us  indicate  certain  properties  of  the  most  simple  linear 
random  functions:  fan  [type]  and  uniform.  .. 

A  fan  [type]  linear  random  function  has  a  common  point  for  all 
realizations  —  the  point  is  called  the  pole.  The  position  of  each 
realization  depends  upon  its  single  random  variable  --the  angular 
coefficient  B. 


-  5 


A  fan  [type]  linear  random  function  nay  be  written  as 

Hit)  '  Hf  +  BCt-tfi) 


(2*13) 


Frequently,  the  abscissa  of  the  pole  tp  is  negative  (Figure  2). 

The  characteristics  of  a  fan  [type]  linear  random  function  are 

W«)  (2.14) 

fty  (*,t')  s  (2.15) 

Assuming  in  formula  (2.15)  t  b  t' .is  the  ratio  between  the 
RMS  deviations 

rnlt)  (2.16) 

Thus  for  fan  [type]  linear  random  functions,  the  RMS  deviations 
depend  linearly  upon  time. 

Normally  a  distributed  fan  [type]  linear  random  function  is 
determined  by  four  numerical  characteristics:^  p,  tp,'’?  t>> 

A  uniform  linear  random  function  is  characterized  by  a  con¬ 
stant  angular  coefficient  for  all  realizations  (Figure  3).  The  ini¬ 
tial  ordinate  Hq  of  a  straight  line  is  the  only  random  quantity 

//(t>  srie  +bt  (2.17) 

Therefore,  one  can  also  represent  a  uniform  linear  random 
function  as  the  random  quantity,  the  center  of  dissipation  of  which 
shifts  in  accordance  to  linear  law  with  the  passage  of  time.  The 
characteristics  of  a  uniform  linear  random  function  have  the  form 

»?„«)  -  »*>•»  (2.18) 

K  (ti')~Dh0  -  ton  ft:  (2.19) 

Thus,  the  constancy  of  its  dispersion  is  the  distinctive 
attribute  of  a  uniform  linear  random  function.  According  to 
formulas  (2.18)  -  (2.19),  the  normally  distributed  uniform  linear 
random  function  is  completely  determined  by  the  three  numerical 
characteristics:  m^  ,  D and  b. 

A  uniform  linear  random  function  possesses  an  important  . 
property  for  practice:  if  one  plots  parallels  to  the  coordinate 
axes  of  a  section  through  the  given  point  k  for  the  expected  value 
of  the  function  (Figure  3),  then:  . 

i  a)  The  coordinants  of  point  k  will  be  the  expected  values 

of  the  random  quantities  which  correspond  to  the  sections 

(2.20) 

b)  The  dispersions  of  the  fandom  quantities  in  the  sections 
will  be  associated  with  the  relationship  , 

D«K:bxDn  (2.21) 

This  property  follows  because  in  the  considered  sections,  the 
random  quantities  and  Tfe  are  association  with  the  linear  relatidn- 


ship 


’b(T 


U) 


(2.22) 
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Inagreeraent  to  (2.22),  when  the  law  of  distribution  is 
normal,  the  law  of  distribution  will  also  be  normal. 

f  The  possibilities  of  the  application  of  linear  random  function 
are  not  limited  by  questions  of  shelf-life  (reliability).  They  can 
be  found  useful  during' the  solution  of  the  series  of  applied  problems. 

,3.  Computation  of  the  law  of  distribution  for  the  hold-up  time 
of  elements  in  proper  condition.  An  element  is  considered  operable 
if  the  value  of  its  determining  parameter  is  greater  than  the  critical 
value  kJ  At  the  fixed  moment  of  time  t^,  .the  probability  g  ^j;  is 
that  the  element  is  operable  (that  is,  that  the  value  of  its  determin¬ 
ing  parameter  Hj  is  greater  than  t*)  )  is  expressed  by  the  formula 

'  {  h  (yi)  j  (3.1) 


where  f^(h  )  --  is  the  probability  density  of  values  of  the  determin¬ 
ing  parameter  at  the  "i"  section.  '  - 

If  the  random  function  is  distributed  normally,  then 

/  -  r_  (n 


2<r„  * 


(3.2) 


Here  m^j  --  is  the  mathematical  expectation  of  the  values 
of  the  determining  parameter  at  the  "i"  section,  and  cr-^  --  is  the 
RMS  deciation  of  values  of  the  determining  parameter  from  the  mean 
value.  ••  y 

One  can  rewrite  formula  (3.2)  as 

£  (<-.)  <3-3> 


where 


"(3.4) 

■  <r»L  .  '  ■  V 

where  £  (nt)  --  is  the  function  of  Laplac.  . 

The  probability  that  the  element  will  go  out  of  order  during 
the  time  (ti,  t^*  »  t^  +  6  t)  will  be  equal  to  a  decrease  during 
this  time  of  the  probability  that  the  element  is  operable 

ft  (6  i)  >,£l  s 

Substituting  the  increment  A  by  the  differential  when 
A  t  0,  we  have  _  ....  .  „  \.L  ; 

f  l  (d  t)  -  -  ii  GiVy  -  cl  /?  UFJ  (3.5) 

where  R«<;  *»  1  -  Gw^  —  is  the  probability  that  the  element  is  not 
operable  at  the  moment  of  time  tj^  if  the  critical  value  of  the 
determining  parameter  is  equal  to  k)  ... 

On  the  other  hand,  the  probability  [of  the  element]  going  out 
of  order  during  the  time  (t^  .ti  +  dt)  is  equal  to 

9:  *  f  I dtj  >t,  (3.6) 

f*/  (t)  --  is  the  density  of  probability  of  hold-up  timeof  the  element 
in  operable  condition  if  its  limit  is  equal  to  o) . 
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Frdm  (3.5)  and  (3.6),  we  have 

fjt)l  •■■■■  . 

'i**i  ,  ci  t f  t ** 


dj«jt)  I 

Jt  1  c*t, 


(3.7) 


Thus,  the  density  of  probability  of  the  hold-up  time  in  oper¬ 
able  condition  fy  (t)  (or,  in  other  words,  the  density  of  probability 
for  the  time  of  a  failure  appearing)  is  the  rate  of  change  of  the 
probability  that  the  element  is  not  operable. 

In  accordance  with  (3.7),  the  probability  of  a  failure  appear¬ 
ing  during  the  time  d  tA  =  t*;  +  1  - , tj  is  equal  to  a  decrease  during 
this  same  time  of  the  probability  that  the  element  is  operable  (or 
to  an  increment  of  the  probability  that  the  element  is  not  operable) 

"  H?i*t  k  (3.8) 

In  accordance  with  (3.7)  -  (3.8),  the  sh«lf-life  of  an 
element  depends  upon  the  nature  of  the  change  of  the  value  Goi;  with 
time,  but  the  value  of  the  probability  GuJl  that  the  element  at  the 
moment  of  time  t^  is  operable  still  does  not  determine  its  shelf -life 
(reliability).  In  other  words',  one  can  speak  only  about  the  shelf- 
life  (reliability)  of  an  element  in  that  case  when  the  segment  of 
time  of  its  existence  (t^,  t±  +  1)  is  .indicated. 

When  t  =  0,  that  is,  at  the  moment  of  issuance  from  the 
factory,  the  element  will  also  be  operable  with  the  certain  proba¬ 
bility  .  At  the  factories,  a  series  of  measures  are  taken  for 

decreasing  the  initial  dispersal  of .  the  determining  parameter's 
values.  Therefore,  the  value  of  probability  G«.^  is  close  to  one. 

In  accordance  to  the  value  Go>i,  ,  one  can  judge  concerning  the 
homogeneity  of  the  product  produced  by  the  factory;  however,  nothing 
can  be  siaid  concerning  its  reliability.  It  is  necessary  to  know 
the  characteristics  of  the  random  function  H(t)  in  order  to  judge 
the  reliability  of  elements  of  a  given  type. 

When  computing  the  density  of  probability  of  the  hold-up 
time  of  an  element  in  operable  Conation,  it  is  necessary  to  solve 
a  quite  unusual  problem  for  the  theory  of  random  functions:  in 
accordance  to  the  characteristics  of  the  random  function  H(t),  one 
finds  the  law  of  distribution  of  a  random  function  at  one  of  the 
horizontal  (that  is,  parallel  to  the  axis  t)  sections.  Since  the 
characteristics  of  the  random  function  H(t)  are  computed  in  accor¬ 
dance  to  the  values'  of  the  random  quantities  H^,  H^,  ...  in  the 
vertical  sections,  then  the  problem  virtually  consists  of  determining 
the  law  of  distribution  of  a  random  quantity  in  one  of  the  possible 
horizontal  sections  according  to  the  values  of  the  numerical  char¬ 
acteristics  of  a  minimum  number  (two)  of  vertical  sections. 

Upon  actual  computation  of  the  density  of  the  probability  of 
hold-up  time  of  elements  in  operable  condition,  one  can  proceed  in 
the  following  manner.  The  entire  range  of  time,  proceeding  from 
the  moment  of  manufacture  of  the  element,  is  divided  into  intervals. 
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For  each  interval  a  tj.  =  tj.  +  i  -  tj.  is  found  the  average  value  of 
the  density  of  probability  of  the  hold-up  time  in  operable  conditions 
f  according  to  the  formula  which  results  from  (3.8) 


‘  j  w  •  i:  iih).  :s  ~  ’'ujl  (3,9) 

•  Xutt  ti  %-i  *1  t  i  ..  tf- 

If  the  random  function  H(t)  is  distributed  normally,  then  one 
can  express  formula  (3.9)  in  accordance  with  (3.3)  and  (3.4)  via 
Laplac's  function 

r  .  <§(«;+,)  -  &(“.)  ' 

ti*.,  -  it  (3.10. 

where  ui  +  i,  u^  are  determined  in  agreement  With  (3.4) .  A  histogram 
is  constructed  in  accordance  to  the  obtained  values  of  f<J«  which 
form  a  continuous  curve 

Thus ,  to  compute  the.  average  value  fit/*  ,  which  corresponds 
to  the  interval  (V  ti*  it  is  necessary  to  know  the  laws  of  distribution 
of  the  element  's  determining' parameter  at  the  beginning  and  end  of  this 
time  interval.  For  the  iiormelly  distributed  random  function  H(t), 
it  is  sufficient  to  know  only  its  characteristics  mjj  (t)  and  (t) 
according  to  which  the  values  mb'  and  cr*  ,  which  correspond  to  the 
beginning  and  end  of  each  of  the  intervals  &  tf,  are  found. 

The  characteristics  of  the  linear  random  function  H(t)  may 
be  determined  from  a  test  according  to  t^e  results  of  the  measured 
values  of  the  determining  parameter  of  similar  elements  at  the  two 
moments  of  time  t^  and  tj_  +  Moreover,  the  relationships  myj  (t) 
andcr’»\(t)  can  be  found  according  to  formulas  (2.2),  (2.4)  and  (2.8) 

-  (2.i2).  ,  ■;  :  .. ...  .. . 

Thus,  the  probability  of  a  hold-up  time  in  operable  condition 
may  be  determined  by  utilization  of  formula  (3.7).  For  the  normally 
distributed  random  function  H(t),  we  have  after  conversion 

f m  /-tVig)  yi"l  6.11) 

According  to  general  formula  (3.11),  one  can  find  the  laws 
of  distribution  of  the  hold-up  time  in  operable  condition  for  various 
individual  types  of  normally  distributed  random  functions;  Having 
substituted  in  formula  (3.11)  the  values  m  t\  (t),^^  (t),  and  uJ  in 
accordance  to  (2.14)  and  (2.16),  and  having  designated  tn ^  /o£  =  c, 
we  will  obtain  after  conversion,  the  formula  of  the  density  of  prob¬ 
ability  of  the  element's  hold-up  time  in  operable  condition  for  the 
normally  distributed  fan  [type]  linear  random  function  H(t) 

■f-  uy.  7  „ 

' '  •:??(“)*  1  i  <3'12> 

where  t^  —  is  the  abscissa  of  the  point  of  intersection  of  the 
straight  lines  m^(t)  and  -  const. 


■fwLi)  = 


For  uniform  linear  random  function,  the  moment  of  time  t„j  will 
be  the  expected  value  (the  average  value)  of  the  element's  hold-up 
time  in  operable  condition  when  the  limit  of  the  element  is  equal  toiO. 

For  a  normally  distributed  linear  random  function  H(t),  making 
a  substitution  of  the  variable  in  formula  (3.11)  according  to  (2.18) 
and  (2.21),  we  will  obtain  after  conversion  the  normal  law  of  dis¬ 
tribution  of  an  element's  hold-up  time  in  operable  condition 


KM)  - 
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Lt'tJS 


x 


(3.13) 


where  is  the  RMS  deviation  of  the  element's  hold-up  time  in 

operable  condition  from  the  average  value  for  the  critical  value 
of  the  determining  parameter  which  is  equal  to  <*>  . 

Thus,  when  the  normally  distributed  linear  random  function 
H(t)  is  uniform,  computation  of  the  density  of  probability  of  the 
hold-up  time  in  operable  condition  results  in  detetminig  t^j  and<r£. 
from  the  very  simple  equations  (2,18)  and  (3.21).  Therefore,  one 
can  raise  the  question  concerning  [whether)  the  substitution  of  the 
available  normally  distributed  linear  random  function  is  equivalent 
to  its  uniform  linear  random  function.  We  will  conditionally  call 
such  a  substitution  the  equalization  of  linear  random  function  H(t). 

Upon  substitution  of  the  linear  random  function  H(t)  of  the 
uniform  linear  random  function  H*(t),  an  error  arises  in  the  cal¬ 
culation  of  the  density  of  probability  of  the  hold-up  time  in 
operable  condition.  It  is  necessary  to  choose  the  parameter  of  the 
univorm  linear  random  function  R*(t)  so  that  this  error  wbuld  be 
minimal.  Overall  consideration  of  the  question  concerning  equalization 
of  a  linear  random  function  goes  beyond  the  confines  of  this  article. 
Let  us  consider  only  one  individual  occurrence. 

When  special  accuracy  is  not  required  or  when  the  random 
function  H(t)  is  dose  to  the  uniform  [function]  H*(t)  (that  is, 
the  dispersion  D  (t)  changes  little),  one  can  assume 


b  *  =  tvtb 


(3.14) 


and  the  RMS  deviation  orv  of  the  element’s  determining  parameter 
is  considered  constant  and  equal  to  its  mean  value  during  the  varia¬ 
tion  of  the  parameter  from>?^  to^^j 


'fij 


jL 

*  z 


fT 


y  -  t  oiis  f" 


(3.15) 


Moreover,  the  numerical  characteristics  tuj  and  of  the 
normal  law  of  distribution  of  the  hold-up  time  in  operabld  condition 
can  be  approximately  determined  by  two  measurements  of  the  parameter 
Y(  in  the  storage  process.  The  value  t,j  is  found  from  the  equation 

uf  -  *  +1>  *  + 


10  - 


(3.16) 


Having  substituted  the  values  b*  and  according  to  (2.8) 

and  (2.9)  in  (3.16),  we  have 


f  -  (-t  »>■  **  '^V/u/  —  X t  j.( J*int  ~t~  )l4/ 

rw  ~  nr  - . - 


(3.17) 


where  ,  nVv*i  are  the  mean  values  of  the  element's  determining 
parameter  at  the  moment  of  time  t^. 

Extrapolating  linearly  the  RMS  deviation  of  the  element's 
determining  parameter  from  its  mean  value,  we  have 


CTk  » 


<71 .  =  <rr  , 
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o  4  j  ~,‘t  i  .  r  ,  . 

•  Substituting  the  value  0»i  wj  and  O'*  ^  in  formula  (3.15)  and 
taking  into  consideration  (2.21),  we  will  obtain 

cT  =•  ^  -  T^t  )  t  <*>  >  4  .  (3.18) 

£  -Z>)  . 

Having  computed  the  values  t  and  3"*  according  to  formulas 
(3.17)  and  (3.18),  one  can  find  the  probability  of  failure  (1  -  s) 
of  the  element  during  storage  time  t  according  to  formula  [3] 

/-  5  =  ^  (3-190 

where  -y  _  y 

(3.20) 

but (Vj  —  is  Laplac's  function. 

The  shelf-life  of  the  element,  in  accordance  to  (3.19), 

will  be  i  jc  /J\ 

5  =  r  +&  M  (3.21) 

Setting  several  values  of  t,  one  can  compute  and  plot  the  shelf-life 
function  s(t),  according  to  which  one  can  judge  concerning  the  possib¬ 
ility  and  the  expediency  of  the  equipment's  future  storage. 

If  in  the  process  of  storage,  there  is  the  opportunity  to 
repeat  measurement  of  the  values  of  the  elements'  parameters,  then 
the  above  stated  method  can  be  applied  repeatedly.  Moreover,  by 
means  of  successive  approximations,  one  can  obtain  still  more 
accurate  computational  results. 
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Figure  1,  One  of  the  methods  for  measurement  of  the  realized  values 


of  a  linear  function. 


Figure  2.  A  fan  [form]  linear  random  function. 
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Figure  3.  A  uniform  linear  random  function. 
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